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ABSTRACT 

We derive accretion rate functions (ARFs) and kinetic luminosity functions (KLF) for jet- 
launching supermassive black holes. The accretion rate as well as the kinetic power of an 
active galaxy is estimated from the radio emission of the jet. For compact low-power jets, we 
use the core radio emission while the jet power of high-power radio-loud quasars is estimated 
using the extended low-frequency emission to avoid beaming effects. We find that at low 
luminosities the ARF derived from the radio emission is in agreement with the measured 
bolometric luminosity function (BLF) of AGN, i.e., all low-luminosity AGN launch strong 
jets. We present a simple model, inspired by the analogy between X-ray binaries and AGN, 
that can reproduce both the measured ARF of jet-emitting sources as well as the BLF. The 
model suggests that the break in power law slope of the BLF is due to the inefficient accretion 
of strongly sub-Eddington sources. As our accretion measure is based on the jet power it also 
allows us to calculate the KLF and therefore the total kinetic power injected by jets into the 
ambient medium. We compare this with the kinetic power output from SNRs and XRBs, and 
determine its cosmological evolution. 
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1 INTRODUCTION 

There is increasing support for the idea that feedback from active 
galac tic nuclei (AGN) plays an important role for galaxy formation 
(e.s., ICowie & Binn ev 1977; Binnev & Tabor 1995i: Isilk & ReesI 
Il998l: IChurazov et al...2001. : .Pi Matteo et al., 2005.) . This feedback 
is also inv oked to explain the M-u relation of supermassive black 
holes fe.g..'Haehnelt et al.ll998l:lKing|2003l : lRobertson et alj2006l : 
iFabian et al. 2006). It is not yet clear whether kinetic or radiative 
feedback is dominant and how efficient each is at a given accretion 
rate. In this paper we will exploit the analogy between X-ray bi- 
naries (XRBs) and AGN to obtain information about the efficiency 
of the different feedback processes and calculate the total power 
available for feedback at a given redshift. 

The central engines of XRBs and AGN seem to be similar 
(e.g.. iMirabel & Rodrigu"e3ll998l : lMeiedl200lh and recently rela- 
tions have been found which sca le spectral and variability proper- 
ties f r om one class to the o ther ( Merloni et al. 2003|; iFalcke et al.l 
|2004 lK6rdingetaDl2006cl : iMcHa rdv et al. 2006). In the nearby 
universe {z < 0.2), there are very few high-luminosity quasars 
like those that exist at high redshifts. These bright high redshift 
quasars are likely to have a strong effect on the evolution of the 
galaxy luminosity function. However, there are several XRBs that 
go through transient phases of very high accretion rates. These 
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"very-high state" objects may be better templates for the bright 
quasars at a redshift of two than any nearby AGN. Thus, we will 
use our knowledge of XRBs and their states to obtain information 
about the kinetic and radiative properties of AGN. 

For XRBs one can observe a full outburst cycle, in which 
the accretion rate increases from very low (10~^ ^5 ^ 



10~^) to near the Eddington limit and then decays again. At 
low accretion rates, the source is generally found in the hard 
state, which is characterized by a hard power law in the X- 
ray spectrum. The hard X-ray emission is usually ac compa- 



nied by radio emission associated with a compac t jet jFendei 



'200 ll) which sometimes can be directly imaged ( [Stirling et al 
2001,, e.g.,). The accretion flow in the hard state is likely to 
be radiatively inefficient (e.g..|Esin. McClintock & Naravan|[T997l : 



iKording. Fender & Migliarj|2006d) . The source can stay in the 
hard state to fairly high luminosities (~ 30% Eddington) until it 
changes its state. During the state transition, it first enters the hard 
intermediate state (IMS), which is characterized by a hard spec- 
tral component and band-limited noise in the power s pectrum. It is 
usual ly accompanied by an increasingly unstable jet lender et al.l 
l2004h . After the hard IMS it enters the soft IMS, which is dom- 
inated by a soft spectral component and power-law noise in the 
power spectrum. Near this transition one often observes a bright 
radio flare; however, after the flare the jet seems to be quenched. 
After the soft IMS, the source is often found in the soft state where 
the X-ray spectrum is dominated by a soft multi-colour black body 
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component. Here the je t is typically quenched by a factor > 5 in 
GHz radio luminosity jpender et al]|l999l : ICorbel et alj|2000l) . In 
the soft state it is generally assumed that the source is efficiently 
accreting, i.e., it has a standard geometrically thin, optically thick 
accretion disk ( Shakura & Sunyaev 1973). On the further way back 
down to low accretion rates, the source stays in the soft state until 
it reache s a critical accret ion rate around ~ 2% of the Edding- 
ton rate iMaccaron3l2003h , where the reverse state transition be- 
gins. This proceeds again via the soft IMS and the hard IMS to 
the hard state. The transition from the hard to the soft state typ- 
ically occurs at a higher luminosity than the transition from the 
soft to the hard state. Due to this hysteresis effect there is no one- 
to-one correspondence of accretion rates to accretion state s. For a 
detailed discussion of state s and their exact definition s see INow^ 
(Il995h:l Belloni et al. I(l2005hj Homan & Bellonil ( [200i) . but see also 
iMcClintock & Remillard booeh for slightly different definitions 
(mainly concerning the intermediate states). 

Both in X-ray binaries and AGN, the unbeame d radio lumi- 



nosity can be used to estimate the accretion rate. Kor ding et al 



2006bh do this for the core radio emission, while iWillott et al 



1999h present a correlation between extended radio luminosity 



and tracers of the accretion rate of luminous double-lobed radio 
sources. Here, we use these methods to translate radio luminos- 
ity functions into accretion rate functions (ARF) of jet-emitting 
sources. 

In order to turn a radiative bolometric luminosity func- 
tion (BLF) into an ARF as well, it is necessary to know the 
radiative efficiency of the accretion flow. There have been 
several suggesti ons that i neffici e nt accretion is visible both 

Jl997h 



in X RBs (e.g., 
(e.g., iRees et aU 



Ichimaru 

1982; Chiaberge et al. 200^, 
of authors have attempted to link both source class e s and 
establish a detailed corresponde nce (see, 
Livio. Pringle & Kinj [2 003; Mac carone. Gallo & Fende 
Merloni, Heinz c& Di Matteo 2003; 'Falcke, Kording & Markof^ 
2004; Jester 2005; Kording, Falcke & Corbel 2006a; Kording et al. 
2006bj). These sources are not only likely to be inefficiently accret- 



and AGN 
and a number 



ing, but the total energ y output of the sou r ces may be dominate d 
by the jet power (e.g.. Fender etai]|2003l ; iKording et al]|2006bh . 
To compare the ARF with the BLF we will explore the effect 
of inefficient accretion flows for low-luminosity o bjects that are 
included in the bolometric luminosity function of ( IHoDkins et al.l 
120071) . Comparing our radio-derived ARF obtained in this way 
with this BLF, we can also determine a radio-loud fraction that 
has an intuitive physical meaning: the ratio of volume densities of 
radio-loud and radio-quiet sources at a given total accretion rate. 

As the jet power is employed in our method for estimat- 
ing the accretion rates from radio data, we can also directly es- 
timate the jets' kinetic power from the radio luminosities (for an 
approach estimating jet pow ers from modeling blazar SEDs, see 
iMaraschi & Tavecchioll200lh . This allows us to construct kinetic 
luminosity functions for jet emitting sources in the local universe 
as well as for higher redshifts. Such kine tic luminosity functions 
have already been constructed bv lHeinz"e t al. (2007). However, we 
use a slightly different methodology and are able to include high- 
luminosity sources (i.e., FR-II radio galaxies). By integrating the 
kinetic luminosity functions over all luminosities we then estimate 
the total power available for kinetic and radiative feedback at a 
given redshift. We also consider constraints on the Eddington ra- 
tio distribution of accreting black holes which reproduce the ob- 
served luminosity functions of all classes of AGN using only the 
black-hole mass function and prescriptions for the radiative and ki- 



netic out puts of accreting black holes as function o f accretion rate 
(compare IVolonteri et al i 20061 ; iMarulIi et al.ll2007h . 



2 ACCRETION RATE AND JET POWER ESTIMATES 
BASED ON RADIO LUMINOSITIES 

Low-luminosity AGN usually have a comp act flat-spectrum radio 
core and relatively weak optical emission fad 1 19991 ; iNagar et al.l 
I2OO5I) . They probably have no "big blue bump", and their total 
energy output of the source may be dominated by the jet. In the 
context of the analogy between XRBs and AGN these objects are 
likely to correspond to hard-state objects, which exist predomi- 
nantly below accretion rates of ~ 2% of the Eddington rate. High- 
luminosity radio sources with significant extended emission are 
typically classified as FR-II radio galaxies or radio-loud quasars. 
These sources do have a big blue bump and are therefore likely 
to be in the analogue of intermediate- state (IMS) sources. The to- 
tal power output of such sources is likely dominated by radiation. 
As the accretion states depend on the Eddington-scaled accretion 
rate (albeit with a hysteresis), a black hole that has a low mass 
but is accreting strongly (e.g., an FR-II radio galaxy) may have 
a lower jet power than a much more massive BH accreting in a 
strongly sub-Eddington regime. Thus, simply referring to the two 
classes of sources as "low" and "high" power sources is mislead- 
ing. We will therefore refer to those sources accreting at a low frac- 
tion of the Eddington rate and showing a compact radio core as 
"LLAGN-like" jet sources. The radiation-dominated strongly ac- 
creting sources with high-power jets will be referred to as "RLQ- 
like" jet sources. This division in two classes of jets is supported by 
similar differences in XRB jets. Hard-state objects have relatively 
stable and compact jets with a flat radio spectrum, while the IMS 
sho ws ejections of highly r elativistic blobs, i.e., an unstable jet. 

iKording et al.l j2006bh have presented a prescription for es- 
timating accretion rates of hard-state XRBs and unbeamed AGN 
from the core radio luminosity. Analysis of the V LBI monitor- 
ing of a sample of flat-spectrum radio sources by ICohen et al.l 
1 I2OO7I the MOJAVE survey) has shown that only high-luminosity 
sources have strongly beamed jets, while low-luminosity sources 
have a lower maximum Lorentz factor, i.e., slower jets. This is 
also implied if the analogy bet ween XRBs and AGN extends to 
jet speeds, as lGallo et al.l ( l2003h found that hard-state jets are un- 
like ly to be stron gly beamed while IMS jets may be faster (but 
see lHeinz& Mer loni 2004). Thus, low-power LLAGN-like AGN 
are probably not affected strongly by beaming and we can use 
the accretion-rate measure from , Kording et al. ( 2006b,) to estimate 
their accretion rates from their core radio power, and hence the ac- 
cretion rate function (ARF) from their luminosity function. We note 
that some FR-I RGs show high apparent velocities, in fact, they are 
thought to be the parent population of the highly beamed BL Lac 
objects. Nevertheless, they usually do not show a strong "big blue 
bump", so they may belong to the LLAGN-like jet sources. How- 
ever, these FR-I RGs with high jet speeds seem to occur mainly at 
high radio luminosities (> 10^'' W Hz"^ Cohen et al. 2007) while 
we are mainly interested at lower radio luminosities. Additionally, 
there do not seem to be any reported LLAGN blazars. While we 
cannot rule out that none of our LLAGN-like sources are beamed, 
the aforementioned arguments suggest that our LLAGN sample 
will unlikely be strongly affected by beaming. 

However, at higher radio powers, it is currently not possible 
to deboost a flat-spectrum core luminosity function reliably and 
convert it into an ARF, since this would require the exact knowl- 
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edge of the distribution of Lorentz factors for a given accretion rate. 
We therefore estimate accretion rates using the core radio emission 
only for low-power LLAGN-like jet sources. For high-power RLQ- 
like jets, we instead use the unbeamed extended low-frequency 
emission of the lobes, whose strength has been shown to correlate 
with the emissio n-line luminosity, and hence accretion rate, by the 
seminal work of lRawlings & SaundersI il99lh . 



2.1 Accretion rates and jet powers from core emission of 
LLAGN 

2.1.1 Accretion rate 

iKording et al.l ( l2006b !) present an estimate of the accretion rate of 
hard state XRBs and unbeamed AGN from the core radio luminos- 
ity: 

12/17 



M « 4 X 10' 



1030 



erg s" 



(1) 



This accretion measure is normalized using XRBs. As we extrap- 
olate this to AGN the uncertainty increases and we may underesti- 
mate the fluxes slightly (see the mass term of the fundamental plane 
in|Kording et al. 2006b). Thus, as a consistency check we will com- 
pare this accretion measure and jet power estimate to other AGN 
estimates (sect. |2.3| below). 



2.1.2 Jet power 

In addition to the accretion rate, one can also estimate the jet's ki- 
netic power form the core radio flux. In fact, the accretion rate is de- 
termined from the core radio power by assuming that the jet power 
is linearly coupled to the accretion rate. The fraction gj of the ac- 
cretion power which is injected into both jets (Pjot ~ V Mc?) 
is nearly a free parameter (qj ^ 1). We now determine this param- 
eter by comparing our accretion measure with a sample of sources 
with jet powers es timated independently of the core power (see also 
iHeinz et al.ll2007l) . 

The X RB Cyg X- 1 has a well-measur ed jet power and core 
radio flux jGallo et alj|2 005: Russell et al] 120071) . We assume a 
distance of 2.1kpc ('Massey et al. 1995) and a typical radio flux 
of 15mly (e.g.. ,Gleissner et al.. ,2004). For the AGN, we use the 
jet powers determined bvlAllenetaiU2006h and core radio fluxes 
from the samples of iGiovannini et ak I l l2005l) ; iNagar et all ( |200^ . 
All of these jet power measurements are based on the power needed 
to inflate an observed "bubble". They are therefore independent of 
our accretion measure based on the radio core flux. 

In Fig. [T] we show the jet powers determined for Cyg X-1 
and the lAllen et al.l jets in comparison to our accretion measure. 
The solid line assumes that the power in one jet is ~ 1/10 of the 
available accretion power. Thus, if the kinetic jet power estimates 
from bubbles are correct, we find gj ~ 0.2 for the entire twin-jet 
system. This is consistent with the rough estimate of the jet power 
given in i Kording et all ( l2006hl) . The total jet power can therefore 
be estimated: 
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Recentlv.'Binn ev. Bibi & Omm a^ ('2007') have reported simula- 
tions which indicate that jet powers derived from the pressure and 
volume of jet bubbles underestimate the true energy input by a fac- 
tor of order 6. In this case the jets would carry at least half of the 
power liberated in the accretion disk. The accretion disk could then 
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Figure 1. Jet power as a function of the accretion rate. We show the power 
injected into one jet as a function of the accretion rate as measured from 
the core radio emission. For Cyg X-1 and FR-I RGs we plot the jet power 
directly estimated from bubbles. For the FR-II RGs of the 3CRR sample we 
plot 0.1r?Mc^ for both methods to estimates the accretion rate. 



radiate only half of the accretion power. IF this is correct, all our 
accretion measure estimates remain applicable, only the jet power 
and any kinetic power estimate presented below has to be multi- 
plied by a factor 6. 



2.2 Accretion rates and jet powers from extended radio 
emission 

2.2.1 Accretion rates 

The bulk of the radio luminosity of powerful radio galaxies is emit- 
ted by extended lobe s with unbeamed, steep-spectrum emission. 
^Rawlings & SaundersI ( Il99ll) have shown that the strength of this 
extended emission correlates with the narrow-line luminosity of 
the optical core. They used the correlation to argue for a common 
mechanism powering both the radio jets and the optical continuum 
and line emission. The reality of this correlation has been confirmed 
by many later studies. In particular. Iwillott et alj ( fl99S ^ presented 
a refined analysis based on complete radio surveys with vastly dif- 
ferent flux limits, ruling out distance effects as origin of the cor- 
relation. They found a correlation between rest-frame radio power 
at 151 MHz and narrow-line luminosity, and concluded that the jet 
power is within one order of magnitude of the radiative power of 
the accretion disc. As the disc luminosity directly measures the ac- 
cretion rate, we can use the correlation to determine an accretion 
rate from 151 MHz radio luminosities for high-power jet sources. 

To allow a direct comparison of the ARF generated from these 
"radio lobe" accretion rates with the ARF derived from the bolo- 
metric luminosity of all quasars {i.e., including the majori ty of 
quasars without strong radio lo bes) by Hopkins et al.l ( 120071) . we 
do not directly use the result of Willott et alj (11999") but renormal- 
ize this relation usin g _B-band lumin osities together with the bolo- 
metric corrections of iHopkins et alj . 

We construct a sample of broad-line quasars with measured 
74 MHz rad io luminosities by c ross-matching quasars from the 
SDSS DR5 (Schneider et al. 2007) with the VLA Low-frequency 
Sky Survey (VLS^l .Cohen et al..2006„.2007) using a matching ra- 
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Figure 2. Nomializing the jet power measures from radio luminosities with 
sources of measured jet powers. On the ordinate (y-axis) we plot the 74 
MHz radio luminosity and on the abscissa (x-axis) we give the absolute i- 
band magnitude. The line represents the linear relation fitted to the data (see 
text). 



dius of 20" (the histogram of radial separations between the closest 
matched sources has a local minimum at this matching radius). We 
use 74 MHz fluxes as this frequency is near the target of 151 MHz, 
especially for the large number of the quasars around « R=: 2, and 
as the VLSS provides an easily accessible deep survey of the full 
northern sky. The resulting matched list has a total of 919 en- 
tri es. We apply the i-band emission-line and if-correction given 
in [Richard s et al. (2006, Table 4) and AT-correct the VLSS data to 
74 MHz rest-frame assuming a spectrum oc 

In Fig.|2]we show the 74 MHz luminosity against the absolute 
i-band magnitude. Following the nearly linear correlation between 
the narrow-line luminosity an d the lo w-fre quency radio luminosity 
found by Will ott et al.l ( 1 19991 see also l^llo tt (20QD), we fit a lin- 
ear dependence between the optical i-band luminosity (which we 
assume is proportional to the ionizing luminosity, and hence to the 
narrow-line luminosity) and the radio luminosity to the data: 



log 1/7 



-QAMi + 16.78 



(3) 



To convert the i-band magnitudes to B-band magnitudes, we 
assume a power-law spectrum with a = —0.5 (see Richards e t al.l 
[2OO6.) which yields B — i = 0.3. For the B-band luminosity we 
find: 



log Lb 



log L7 



W 



Hz sr 



+ 18.6. 



(4) 



Throughout this paper we mainly use cgs units. However, as most 
radio luminosity functions are given in WHz^^ sr^^ we provide 
the conversion formulae from radio luminosity to ac cretion rate us- 
ing the se units on the right hand side of the equation. iHopkins et al.l 
use a non-constant bolometric correction of the B-band 
flux. However, the bolometric correction deviates from a constant 
mainly at lower luminosities where we will not use the relation to 
obtain accretion rates. As we would like to obtain a simple relation 
between the low frequency radio luminosity and the bolometric lu- 
minosity we will use a constant bolometric correction of 10, i.e., 
Lboi ~ WLb- W ithin our uncertainties , this is in agreement with 
the value used by iHopkins^ 'Zl ( I2OO7I) . If we again assume that 
radio emission has a spectrum oc v~'''^, we can translate the 
measured relation to: 



log Lb 



s J 



log L15 



W 

Hz sr 



19.8. 



(5) 



To obtain the accretion rates we assume that the sources are accret 
ing efficiently with a constant efficiency of r; — 0.1: 

W 



logM 



[f] 



log L151 



Hz sr 



0.15 



(6) 



The correlation between optical narrow-line luminosity and 
151 MHz radio luminosity jWillott et al.] Il999l) has only been 
tested with a sample in the luminosity range of 10^*'^ ^ 
Li5iMHz/(WHz"^sr"^) < 10^*'*. Therefore, we win only use 
it for extended radio emission brighter than lO^**'^ WHz"\ 



2.2.2 Jet power 

If we assume that the ratio of jet power to accretion rate gj is simi- 
lar for hard-state and IMS objects, gj ~ 0.2, we can use the accre- 
tion rate determined from the extended low-frequency radio emis- 
sion also as a measure of the jet power. 1 Willott et aL (,1999r) and 
lHardcastleetaLl ( l2007l) report that 



rjot 



[?1 



3 X lO'^Lisi 



W 



Hz sr 



r 



3/2 



(7) 



where / parameterizes our uncertainty of the Jet power compared 
to the minimum energy needed to account for the synchrotron 
emission from the lobes. iBlundell & Rawlingsl Jioool) suggest that 
/ « 10 is applicable to FR-II RGs. We find: 



logPjc 



[?1 



18.7 + log Li5i 



W 



Hz sr 



(8) 



Our normalization of the L151 -accretion rate correlation (eq.|5j to- 
gether with qj — 0.2 gives a n ormalization consta nt of 19.1 com- 
pared to 18.7 as estimated from lWillott etaP i ll999l) . The difference 
corresponds to a factor 2.5. This is well within the uncertainties of 
our accretion rate and jet power estimates, but may indicate that our 
normalization is slightly too high. 

For the rest of this paper we will assume that the coupling 
constant of the jet is qj = 0.2 and use this to estimate the jet 
power from the accretion rate for both low-luminosity and high- 
luminosity objects. Thus, our jet power measure from extended ra- 
dio emission is: 



logPjc 



[erg 



19.1 + log Li5i 



W 



Hz sr 



(9) 



2.3 Comparison of jet and accretion rate measures from both 
methods 

In the preceding subsections we have presented two different accre- 
tion rate and jet power measures based on radio luminosities. As we 
will use these to obtain accretion rate functions, it is important that 
the estimates are consistent with each other. 

To compare the accretion rate measures, we need a sample 
which has measured values both for the extended low-frequency 
radio flux and for the unresolved core flux, so that we can compute 
and compare both accretion rate measures for the same sample. The 
FR-II radio galaxy subsample of the 3CRR catalogue fulfills these 
criteria. We take 151 MHz radio fluxes from Laing et al. (1983) and 
core radio fluxes from the compilation by M. Hardcastlfl Since the 
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Figure 3. Beaming Statistics of the core radio measure compared to tlie 
accretion measure based on the extended emission: On the x-axis we show 
log Pjct, core /Jjot, lobe while the y-axis denotes the number of sources. 
The filled circles represent the measured histogram including upper limits 
and the triangles show only those sources with ineasured core fluxes. The 
lines describe a beaming model assuming a Gaussian distribution of Lorentz 
factors around F = 4, 6, 8 with cr = 1. 

core fluxes of this sample are probably affected by beaming, we can 
at the same time assess the impact of beaming. 

In Fig. [T]we show the 3CR R FR-II sample together with Cyg 
X-1 and the lAlIen et al.l ^20061) sample of sources with have jet 
powers inferred from X-ray bubbles. Most of the 3CRR sample 
lies above the expected line, i.e., the jet power estimated from the 
extended emission exceeds that estimated from core radio fluxes. 
While the core flux is affected by beaming, the extended emission 
is not. For a randomly orientated sample of sources, like the low- 
frequency selected 3CRR sample, the majority of sources will be 
(ieboosted, since they are beamed into a narrow cone (for F > 1.7). 
This explains the lower average jet power estimated from the core. 

An additional problem that may alter the relative distribution 
of the two accretion power estimates is that the extended emission 
of the lobes can only yield a measure of the jet power averaged 
over the lifetime of the lobes (several million years). In contrast, 
the core radio fluxes measures the instantaneous jet power. Thus, if 
a source changes from having a strong jet to having a quenched jet, 
but the relativistic particles in the lobes keep radiating for at least a 
synchrotron lifetime, this will reduce only the core power measure. 
The reverse effect, i.e., a starting jet, can also happen. However, as 
the 3CRR sample is selected using low-frequency radio emission 
those starting sources will not be in the sample. 

To verify that the ratios between the core jet power estimate 
and the lobe jet power estimates are consistent, we show a his- 
togram of the measured ratios in Fig. [3] Some of the sources in 
the 3CRR catalogue have no measured core radio fluxes. We there- 
fore present two histograms: In the first, we only include detected 
sources. In the second we include all non-detected sources at their 
detection limit. In addition to the measured distribution we also 
present the distribution expected from a randomly orientated sam- 
ple of relativistic jets. Here, we assume that the observed lumi- 
nosity is relativistically boosted as L ~ D^'^ Lq, where Lq is the 
luminosity in the rest-frame of the jet and D is the Doppler factor 
D — 1/(7(1 — f3cos9)]. We have assumed that the Lorentz fac- 
tors are Gauss-distributed around a Lorentz factor Fo with a = 1. 
We show the model for three different mean Lorentz factors: F — 



Figure 4. Different models for the cosmological evolution of the bolometric 
luminosity function extrapolated to 2: = 0. The points represent the mea- 
sured LLAGN radio luminosity function translated to accretion rates. To 
allow for an easy comparison with the bolometric luminosity function we 
show the accretion rates as an equivalent bolometric luminosity (0.1 Mc^). 
To present the uncertainties due to the model of the cosmological evolution 
we show four different models: Full evolution, PLE, PDE and LDDE. Only 
the latter model is used in the remainder of this paper 

4, 6, 8. While the model assuming F = 4 over-predicts the number 
of boosted objects (Piote/ Pcore 1? 1) it under-predicts the num- 
ber of strongly deboosted sources. The model assuming F = 8 
strongly over-predicts the number of deboosted sources. The data 
is best represented by the model with F = 6. This is consistent 
with the average Doppler factors found in BL Lac objects (D ^ 7; 
iGhisellini et al][r993.) . The deviations can be attributed to the sim- 
plicity of the model, as well as to short-term changes in the accre- 
tion rate which only show in the accretion measure based on the 
core fluxes. We conclude that both accretion rate measures yield 
results which are in agreement with each other. 



3 RESULTS 

We have constructed accretion and jet power measures from ra- 
dio luminosities. We now use these to translate radio luminos- 
ity functions into accretion rate functions and compare those to 
the measured bolometric luminosity function of quasars (BLF) 
from Hopkins et al. (2007). For low -luminosity objects we use the 
local radio luminosity function of iFilho et al.l 12006). The high- 
luminosity end is obtained from the low-frequency radio luminosity 
function of Iwillott et al] feOOlh . 

3.1 Local accretion rate functions 

In Fig. |4] we show the local ARF of LLAGN-like jet sources ob- 
tained from our accretion rate measure together with several mod- 
els of the BLF in the local universe (extrapolated to z = 0, 
iHopkins et ai]|2007h . To compare the ARF with the BLF, we de- 
fine an equivalent bolometric luminosity L^'^ = O.lMc^, i.e., 
we plot the ARF as if the sources were accreting efficiently, with an 
efficiency of converting rest-mass energy into radiation of r; = 0.1. 
If the majority of the sources in the BLF are indeed accreting effi- 
ciently, we can compare the ARF and the BLF directly in this way. 
This assumption is likely to be true for the high-luminosity quasars. 
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Figure 5. Bolometric luminosity function at z = ; we also show the 
measured ARF of "low"-power jet sources (LLAGN-sample) and RLQ-like 
jet sources (e.g., radio-loud quasars). The latter ARF is obtained from a 
151 MHz luminosity function together with eq. ^6). The fit to the "low"- 
power jets is the model of the BLF modified to include an exponential cut- 
off as seen in the black hole mass function - see text. 



We will discuss the effects of inefficient accretion on the shape of 
the BLF in section [43] 

Of the models given in ' Hopkins et al] j2007b . we show the 
pure luminosity evolution (PLE) model, the pure density evolution 
(PDE) model, the luminosity dependent density evolution (LDDE) 
model and the full evolution model. The latter is basically a broken 
power law model where all parameters (both slopes and the break 
location) can have a different evolution (for the exact definitions of 
the different models see the original paper). We show the different 
models mainly to demonstrate the uncertainties of the modeling and 
will adopt the LDDE model for the rest of this paper as it currently 
seems to be the most widely used model and the 151 MHz luminos- 
ity function is fitted with an LDDE type model. The BLF has data 
reaching to 10^^'^ erg/s, below that the models are extrapolated. 

Our ARF obtained from the local radio luminosity function 
is in agreement with the BLF, within the uncertainties of the cos- 
mological evolution. Thus, at low accretion rates there are similar 
numbe rs of radio jet-e mitting sources and sources appearing in the 
BLF of iHopkins etZl 

Our accretion measure using core radio fluxes is strongly af- 
fected by relativistic beaming. However, no AGN below a radio 
luminosity of ~ 5 x 10^^ W/Hz shows high apparent velocities 
dCohen et alj|2007[) . This radio luminosity corresponds to a bolo- 
metric luminosity of « 2 x 10*"* erg/s. This supports the idea that 
the low-luminosity radio LF is not strongly affected by relativistic 
beaming. We mentioned that we expect LLAGN-like jets to exist 
predominantly below ~ 2% of the Eddington rate. For a 10** M0 
BH this corresponds to a bolometric luminosity of 



4 X lO'*'' erg 



s~ . Thus, LLAGN-like jet are likely to be slower than RLQ-like 
jets. 

The majority of the sources in the low-luminosity radio LF 
would not be classified as Quasars as they are faint in the optical 
band a nd/or do not show br oad emission lines. Even though the 
BLF o f 'Hopki ns etallbOOTh is called a "Quasar LF", it also con- 
tains lower-luminosity AGN in addition to Quasars. This partly ex- 
plains why the ARF of LLAGN-like jets extends the BLF smoothly 
to lower accretion rates. This will be further discussed in Sect. 14. II 

In Fig. [5] we show the measured ARF of the LLAGN-like jet 



sources together with the BLF and the ARF deduced from the 151 
MHz radio luminosity function. The latter ARF describes the abun- 
dance of radio-loud Quasars with lobes, i.e., radio-quiet quasars 
will not be picked up. While jets of LLAGN-like sources do not 
seem to be highly relativistic, those of RLQ can be highly relativis- 
tic. We will refer to the ARF deduced from the 151 MHz LF also as 
the luminosity function of "RLQ-like" jets, and we will use the term 
"strongly jet emitting sources" to refer to all jet-emitting sources, 
whether the jet is RLQ-like or LLAGN-like. Sources without strong 
jets will be referred to as radio-quiet. 

While the LLAGN-like jet ARF is roughly in agreement with 
the BLF, the RLQ-like jet ARF lies significantly below it. How- 
ever, continuity demands that the overall ARF of all jet-emitting 
sources connects the RLQ-like jet LF smoothly to the LLAGN- 
like jet ARF. The exact shape of the ARF between the observed 
domains of LLAGN-like jets and RLQ-like jets is unknown. We at- 
tribute the LLAGN-like jet sources to the hard state, which exists 
predominantly below an accretion rate of ~ 2% Eddington. For a 
given accretion rate only black holes above a corresponding mini- 
mum mass will have an Eddington ratio below ~ 2% and always 
have LLAGN-like jets. Above that Eddington ratio it is likely that 
the majority of AGN will either show RLQ-like jets or be radio- 
quiet. Therefore, we can obtain constraints on the ARF of LLAGN- 
like jet sources from the mass function of supermassive black holes 
in th e local univers e. 

IShankar et alj j2004l) report that the mass function can be de- 
scribed by a power law with an exponential cutoff: 



$bh(M) = 7.7 X 10" 



M 



exp 



M 



,(10) 



with M, = 6.4 X lO''A/0. The BH mass function describes how 
many BHs can harbour LLAGN-like jets for any given accretion 
rate. The measured LLAGN-like jet ARF can be well described by 
a power law - similar to the BH mass function at low BH masses 
(albeit with a different power law index). As the BH mass function 
has an exponential cut-off the power law describing the LLAGN- 
like jet ARF has to cut-off at some point as well. If the majority 
of AGN with an Eddington ratio larger than ^ 2% does not show 
LLAGN-like jets, the BH mass function suggests that the cut-off 
of the ARF will be the same exponential translated to a luminosity 
corresponding to ~ 2% Eddington for BHs with mass Af* (this lu- 
minosity is ~ lO'**'^ erg/s). As the LLAGN-like jet ARF is roughly 
in agreement with the BLF, we therefore describe the LLAGN-like 
jet ARF with the analytic model of the BLF with an additional ex- 
ponential cutoff and a constant logarithmic offset. This model with 
an offset of 0.5 dex is shown in Fig.[5] The total ARF of jet-emitting 
sources is the sum of both radio ARFs derived from RLQ-like jets 
and LLAGN-like jets, respectively. 

We can now revisit our assumption that the low-luminosity 
AGN have slow jets. If this assumption was wrong, we would have 
to correct the luminosity function for beaming. However, we are 
mainly interested in the low-luminosity end of the radio luminosity 
function. At the lowest luminosities the beaming cor rection woul d 
in fact increase the total number of sources (see e.g., lListeill2003b . 
However, the ARF already gives a density of ~ 10~^ BHs per 
Mpc'' and log M for accretion rates below an equivalent bolometric 
luminosity of ~ 10'*'^ erg s^^. If the beaming correction increased 
this value significantly (e.g., a factor 10) this would suggest that 
there were more low-luminosity AGN than BHs in the local uni- 
verse. Thus, even if there is a small beamed contribution to the 
luminosity of LLAGN, the beaming correction cannot be strong. 
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Figure 6. Luminosity function at differe nt redshifts: In the le ft panel we show the evolution of the BLF at four different redshifts (0,0.5,1,2). This is a 
reproduction of the LDDE model given in lHopkins et al.l )2007l) to allow for a direct comparison with t he ARF of RLQ-lik e jets. In the right panel we show 
the evolution of the ARF of RLQ-like jet sources. The ARF shown is based on the 15 1 MHz radio LF of lwillott j200lh . 
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Figure 7. Fraction of radio-loud objects (i.e., voluine ratio of radio-loud 
objects to the BLF) as function of accretion rate for various redshifts 
(0,0.5,1,2). At low luminosities all sources are likely to be radio-loud. No 
evolution is visible, as we assume that the radio-loud LLAGN extend the 
BLF at low luminosities. The evolution at high luminosities is low for inod- 
erate redshifts (z < 2). Only at intermediate luininosities there seems to be 
a small region showing some evolution. However, this is likely an artifact 
of the choice of our models. Our accretion rate measures work best at the 
lowest and highest luminosities. 



3.2 Redshift evolution of accretion rate functions 

In Fig. [6] we show the redshift evolution of the BLF (left-hand 
panel) and the evolution of the ARF of RLQ-like jet objects (right- 
hand panel). As the LLAGN-like jet ARF at redshift zero is roughly 
in agreement with the BLF, we assume that this is also the case for 
other redshifts. In this case, its cosmological evolution is identical 
to that of the low-luminosity end of the BLF and we do not show it 
separately. 

In Fig.|7]we show the fraction of strongly jet emitting objects 
for a given accretion rate per volume. The number of strong jet 
sources is the sum of the number of the LLAGN-like and the RLQ- 
like jets. At low accretion rates, nearly all sources are radio-loud. 
The radio ARF is in agreement with the low-luminosity end of the 



BLF, or lies slightly above it. At the high-luminosity end, the num- 
ber of RLQ-like jet sources of a given accretion rate is significantly 
lower than the number of quasars (~ 1 dex). The exact fraction 
of jet emitting sources is strongly dependent on the normalization 
of the accretion rate measure based on the extended radio emis- 
sion. The slope of the BLF is ~ —2.5. Thus, all uncertainties in 
the estimation of the accretion rate enter at least quadratically in 
the radio-loud fraction. We found in Sect. I2.2.2| that our jet po wer 
estimates are slightly higher than those of IWillott et alj ( 119991) . If 
we overestimate the jet power by 0.2 dex, the true fraction of radio- 
loud objects at high luminosities is less than 3%, well in agreement 
with the values typically found. 



3.3 Jet kinetic luminosity functions 

As our accretion measures from radio fluxes are based on the ra- 
dio emission from the jet, they can also be used as a measure of 
the jet power (see Fig. Thus, we can compute the kinetic lu- 
minosity function (KLF) for jets, i.e., measure how much power 
is injected into the interstellar or intergalactic medium (ISM/IGM) 
by jets from sources with a given accretion rate. In sect. 12.2.21 we 
found that roughly a fraction gj = 0.2 of the accretion power is 
injected into the jet for all jet sources. 

In radio-quiet quasars, usually a weak compact core is found. 
This suggests that there is still an active jet, albeit one of lower 
power. A typical quenc hing factor for the radio emission of radio- 
quiet jets is ~ 100 (e.g.. lKellermann e t al. 1989). If we assume that 
the radio emission is quenched by a factor 100, the jet powers in 
i (RQQ) are reduced by ^ 
■ iKording et aLll2006bl) . 
effect of the RQQ in the KLF to that of RLQs and LLAGN-like jet 
sources. As the fraction of radio-loud sources is small at high lumi- 
nosities (a few percent, Fig. IT}, the number of radio-quiet sources 
is approximately the number of all quasars. 

The resulting jet KLF is shown in the left-hand panel of Fig.[8] 
The LLAGN-like jet sources inject significantly more kinetic en- 
ergy into the ISM than both radio-loud and radio-quiet quasars. 
The right-hand panel shows that this is true at all redshifts. This 
can be seen even better in Fig|9l where we plot the total jet power 
integrated over all accretion rates as a function of redshift. The 



in radio-quiet sources (RQQ) are reduced by a factor 30 (assuming 
Pj;,t, e.g.. iKording et alj|2006bl) . We can compare the 
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Figure 8. The kinetic luminosity function of jets sliows tlie kinetic power injected into tlie ISM/IGM per Mpc'' and logio Pjct as a function of jet power Pjot ■ 
The left-hand panel shows the different components of the kinetic luminosity function at redshift zero: LLAGN-like jet sources, RLQ-like sources as well as 
radio-quiet sources. The right-hand panel shows the evolution of the kinetic luminosity function of jet sources (RLQ-like and LLAGN-like jets, solid lines) 
and radio-quiet sources (dashed lines). 
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Figure 9. Total kinetic power injected into the ISM at a given redshift by 
AGN. We show the integrated kinetic energy injected into the IGM by 
LLAGN, by radio loud quasars and by radio quiet quasars as well as the 
power available for feedback during the quasar phase (e.g., as radiative 
feedback). The points denote the kinetic energy injected by supernovae. 
Additionally, we show the total bolometric luminosity, i.e. the total power 
liberated through accretion. 



LLAGN-like jet sources dominate the total kinetic jet power out- 
put by an order of magnitude. The contribution of the RLQ-like 
jet sources is just comparable to those of radio quiet sources. RLQ 
create jets that are roughly two orders of magnitude stronger than 
those of the RQQ, but only a few percent of quasars are radio-loud 
(RLQ-like jet sources). Both factors average out and we find that 
the contribution of the sample of RQQ and RLQ to the kinetic out- 
put is roughly similar. 

In addition to the kinetic power injected by AGN, we show 
the energy available for "quasar-mode" feedback in Fig.|9l assum- 
ing that 5 % of the quasar's b olometric luminosity i s available to 
heat the surrounding gas (see iDi Matteo et aL 1 12001 e.g. ; we dis- 
cuss the possible nature of the coupling between gas and radiation 
in sect. l4.2t . Even including quasar- mode feedback, the population 



of low-luminosity AGN with their LLAGN-like jets dominates the 
total power available for feedback at low redshifts. As the cosmo- 
logical evolution of the high-luminosity quasars is stronger than for 
the low-luminosity end, radiative feedback increases in importance 
until it reaches the same order of magnitude as kinetic feedback at 
a redshift close to 2. We note that it is not yet possible to observe 
an LLAGN ARF at 2: = 1-2. The only hints to the evolution of the 
low-luminosity radio LF come from the fact that the LLAGN ARF 
smoothly extends the BLF to lower luminosities, which suggests 
that one can use evolution of the low-luminosity end of the BLF to 
estimate the evolution of all LLAGN. 

It is now widely accepted that the kinetic energy input from 
supernovae has a profound impac t on the formation and evolution 
of galaxies (see lEfstathioull2000l e.g., and references therein). A 
typical supernova injects around 10^^ erg of kinetic energy into the 
ambient medium (e.g.. Chevalier 1977; Korpi et al. 1999). Assum- 
ing that both core-collapse and type la supernovae create a kinetic 
energy of 10^^ erg, we can estimate the kinetic power injected into 
the interstellar medium by supernovae, assuming the high-redshift 
supernova rates estimated from the "Great Observatories Origins 
Deep Survey" by Dahlen et al, (2004), and compare it to the power 
created by jets of supermassive black holes. In Fig.|9|we present our 
power estimate of jets, radiative processes and supernovae. While 
the contributions of high-power jets (i.e., radio-loud quasars) and 
radiative processes to the heating of the ambient medium are signif- 
icantly below the supernova heating rate, LLAGN-like jet sources 
inject a similar amount of power into the interstellar/intergalactic 
medium as supernovae. Also the jets of stellar accreting black holes 
and neutr on stars contribute to the kinetic power, albeit at lower 
rates (e.g.. lFender et al.ll2005h . 

Of course, the impact of jets on the galaxy formation process 
depends critically on the relative magnitude of the jet power, the du- 
ration of each jetted phase, the collimation and absolute jet power 
(see section |42t , and the binding energy and cooling time of gas in 
the galaxy. Therefore, it would be highly desirable to derive Fig.|9| 
separately for different galaxy types, luminosities and black hole 
masses. This requires knowledge of the radio and bolometric lumi- 
nosity functions separately for each subcategory, as well as reach- 



© 2007 RAS, MNRAS 00Q.[TlfT2l 



Accretion rate and kinetic luminosity functions 9 



ing lower luminosities at higher redshifts than is presently possible, 
and is beyond the scope of this work. 



4 DISCUSSION 

4.1 The break in the BLF and the missing LLAGN 

In Sect. 13.11 we showed that the low-luminosity radio ARF is 
roughly in agreement with the AGN BLF. However, supermassive 
black holes at low accretion rate are likely to be radiatively ineffi- 
cient and would not be classified as a quasar. Assuming the usual 
high radiative efficiency of quasars will then underestimate the ac- 
cretion rate. At the lowest luminosities, this suggests that the radio 
ARF and the BLF should increasingly deviate from each other with 
decreasing accretion rates, with the radio ARF exceeding the AGN 
BLF significantly. This is not what we find in the ARFs we have de- 
rived (Figures|4]and|5}. Here, we explain the observed break in the 
BLF as due to the onset of inefficient accretion and the "missing" 
radio LLAGN as due to the fact that there is only a finite number of 
supermassive BHs. 

We have several constraints for modeling the accretion rate 
function of supermassive black hole s. Fir st, the lo cal b lack-hole 
mass function has been estimated bv lShankar et al.l ^2004) . We can 
therefore write the accretion rate function as: 



$(A/) = / $BH(Ai")$(M,M)dlogM, 



(11) 



where "1>(M) is the measured accretion rate function in units of 
Mpc~"^ dlogjQ . Secondly, as the number of AGN integrated 
over all accretion rates down to has to be equal to the number of 
all supermassive black holes, it follows: 



/ 



$(M,M)dlogM. 



(12) 



The distribution l>(M, M) describes the probability to find a black 
hole of mass M with an accretion rate M. In principle that function 
can be an arbitrary function and there is no unique inversion that 
obtains $(M, M) from eq.[TT]with only the knowledge of $(M) 
and $bh(M). As a first approach, let us assume that the distri- 
bution depends only on m, the accretion rate in Eddington units. 
This assumption is equivalent with a distribution M) that 

separates into a function depending only on M and another one 
depending only on M: $(M, M) = $i(M)$2(M). In astrophys- 
ical terms, this assumption means that the Eddington luminosity, 
i.e., the black-hole mass, sets the total power that is available for 
feedback from a given black-hole mass, independently of the host 
galaxy properties. 

It is generally assumed that strongly accreting objects are ra- 
diatively efficient. As the BLF is well-described by a power law 
at high luminosities (above the break), it is likely that the distribu- 
tion $(m) is a power law at high Eddington ratios. However, as the 
power-law index of the luminosity function is steep (~ —2), it has 
to have a cut-off towards low luminosities so that its integral does 
not diverge. As there are known strongly sub-Eddington sources 
(e.g., Sgr A*), a hard cut-off would not be appropriate, and we use 
a broken power law of the form: 



(rh/rhbr. 



for rhbrcak ^ m 
for m ^ rhbrcak • 



(13) 



As the function has to be integrable with respect to d log rh we re- 
quire 7i ^ 0. At high accretion rates, we truncate the power law at 
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Figure 10. Measured ARF and BLF together with our model for the LFs. 
While the model for the ARF assumes that the accretion rate is directly 
measured we assume inefficient accretion below a critical accretion rate for 
the simulated BLF. We can reproduce both, the ARF and the BLF with our 
simple model. 



10 Eddington rates. The distribution will be normalized according 
to eq. {V2[ . 

While the ARF obtained from the radio estimates should di- 
rectly measure $(M), the bolometric luminosity depends linearly 
on the accretion rate only for the high-accretion rate sources. For 
strongly sub-Eddington sources we have to translate the accretion 
rate to a bolometric luminosity. The analogy of AGN to XRBs sug- 
gests that the bolometric luminosity of an accreting BH can be de- 
scribed as: 



Lhoi ~ 



L^dd rn 

L-Edd Tl, 



for m ^ TTlcrit 

for rh < merit 



(14) 



where rh is the accretion rate in Eddington units and merit ~ 0.05. 
Note that the rhbrcak used in the previous paragraph has no direct 
relation to rhcrit- The exact value of rhcrit seems to be different 
for the rise of each XRB outburst (0.03 ^ rhcrit ^ 0.3). On the 
decline it seems to be more stable (rhcrit ~ 0.02). For our simple 
model we will use a constant critical accretion rate of rhcrit ~ 0.05. 

We can now combine the measured supermassive black hole 
function (eq.llOt with our assumptions about the accretion rate dis- 
tribution (eq.|13l) to obtain an accretion rate function, c[>(A/). This 
function can be compared directly to the ARFs obtained from the 
radio emission. To obtain a BLF we include the description given in 
eq.[T4]to translate from accretion rates to luminosities for sources of 
a given mass. The free parameters are 71,2 and the break (rhbrcak) 
in the power-law distribution of rh, and to some extent the critical 
accretion rate rhcrit . The latter is constrained from the observations 
of XRBs and the finding that this v alue seems to be similar in AGN 
( |jesteJl2005l : lK6rding et alj2006bh . 

In this paper we only try to explain the small difference be- 
tween the measured radio ARF and the BLF at low luminosities. 
Thus, we will not explicitly find a best-fit model to the data, but 
only compare a plausible model to the ARF and the BLF. 

The high-luminosity slope of the ARF and the BLF is mainly 
affected by the power-law index 72 at high accretion rates. The 
position of the break in the accretion rate distribution (rhbrcak) 
strongly affects the normalization of the BLF. If the break is too 
low, most BHs have a very low rh and there are not enough bright 
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quasars. Thus, both 72 and rfibrcak can be estimated from the high- 
luminosity end of the BLF above the break. The low-luminosity 
index 71 has a rather minor effect on both the ARF and the BLF 
as long as 71 > 0. In Fig. [TO] we show the local radio ARF and 
the BLF together with one realization of our model (72 = —1.9, 
n^b^oak = 2x 10~^, and 71 = 1). Both the BLF and the radio ARF 
are reproduced. 

The features in the BLF and ARF can be understood as fol- 
lows. The break in the BLF is created by radiative inefficiency at 
low accretion rates (this effect of radiative inefficiency on the lu- 
minosity function of AGN was exploited bv |jestej200g to look for 
inefficiency in AGN). The flattening of the ARF at lower luminos- 
ity is due to the finite number of supermassive black holes - if the 
ARF did not flatten, there would need to be more low-luminosity 
AGN, not leaving enough black holes to account for the observed 
high-luminosity (quasar) population. 

With the simple assumption that the distribution function of 
accretion rates $(m) depends only on the accretion rate in Edding- 
ton units and can be described with a broken power law we seems 
to be able to reproduce the observed LFs. However, there are sev- 
eral sugges tions that the distr ibution does depend on the black hole 
mass (e.g.. lAIlen et al. 1 1200^ . If we relax our assumption that the 
distribution function $(M, M) separates into distribution for the 
accretion rate and one for the black hole mass, we have enormous 
freedom to choose the distribution. We tried several functions, e.g., 
we explored if the break of the broken power law in the m distribu- 
tion can depend on the mass, but were not able to improve the fits 
sign ificantly above t he simp le model we presented. 

iHopkins et aT I ( I2OO6 j) model the quasar luminosity function 
by considering the birth rate of quasars with a given peak luminos- 
ity and the subsequent evolution of the quasar light curve in their 
merger-driven model. They explain the break in the LF as a char- 
acterist ic luminosity at which the birth rate peaks. iHopkins et al.l 
( l2006hl) examine the evolution of the light curves in more detail 
and show that quasar light curves vary systematically with peak lu- 
minosity, i.e., an individual black hole's accretion rate evolution de- 
pends on the peak accretion rate. In the framework of these models, 
$(m) will vary with redshift depending on the age distribution and 
lightcurves of currently active black holes, and possibly in different 
ways from the variation we would require. A detailed comparison 
between our and the Hopkins et al. models woudld require knowl- 
edge of the high-redshift black-hole mass and LLAGN luminosity 
functions, which are however not yet available. 

We note that our model suggests that the weakest LLAGN do 
not only have a small central BH but should also be fed a t a strongly 
sub-Eddington (and even sub-Bondi) accretio n rate ( I Yuan et al.l 
I2OO3I : ICuadra et alj|2006l ; iMarrone et al.ll2007l. see, e.g., the sim- 
ulations and model-fitting for Sgr A* [). 



4.2 Implications for feedback 

4.2.1 Available feedback power as function of accretion rate 

In sect. |3.3| we have confirmed earlier findings that AGN at lower 
luminosities (e.g., LINERs and low-luminosity FR-I RGs) are 
likely to dominate the total kinetic power output in the local uni- 
verse with their LLAGN-like jets ( iBest et al.l |2006| ; iHeinz et all 
Their total kinetic power is comparable to the power injected 
into the ISM by supernovae. The population of jets in radio-loud 
and quiet quasars carry more than an order of magnitude less power 
than the jets of low-luminosity AGN; in fact their total power is not 
much greater than that injected by jets of XRBs. At higher red- 



shifts the power injected into the IGM by quasars increases more 
strongly than the power injected by the lower-luminosity objects. 
However, the low-luminosity objects still seem to dominate. If a 
fixed fraction of the luminosity of quasars is available for feedback 
even in the absence of strong jets (e.g. by direct radiative heating 
or via radiation pressure-driven winds), feedback from luminous 
quasars also increases with redshift. For a feedback fraction of 5%, 
this "quasar mode" would make a significant contribution to the to- 
tal power available for feedback at redshifts close to 2, and even 
dominate the total feedback for fractions larger than 5%. This is 
certainly only true when summing over the entire population of all 
supermassive black holes; a different scenario is likely to be found 
in the center of clusters and in individual galaxies. Again, we stress 
that radio and bolometric luminosity extending to lower luminosi- 
ties, higher redshifts, and separated by host galaxy type, are needed 
to assess the impact of feedback separately for galaxies of different 
galaxy type, luminosity, or black hole mass. 

Nevertheless, from the models used to explain the measured 
radio accretion rate function, kinetic power function and the bolo- 
metric luminosity, we can deduce feedback properties of individ- 
ual objects. The analogy of AGN with XRBs and our explana- 
tion of the break in the quasar LF in the preceding section suggest 
that all accreting black holes below a critical accretion rate around 
3 X 10~^ of the Eddington rate are accreting inefficiently and in 
a jet-dominated mode. Above that critical accretion rate, a source 
can either be in the hard or hard intermediate state and continue to 
be strongly jet emitting, or the jet power is quenched by a factor 
~ 30. With the efficiencies discussed above, we obtain the follow- 
ing approximations for the kinetic power available for feedback, as 
function of the fundamental parameters m and L-^dd, i-c, Mbh'. 

{0.2m Lsdd if rh ^ riicrit 

0.2 epR-ii mLEdd if m > merit, with jet (15) 
0.007mLEdd if rh > riicrit, no jet 

The "with jet" and "no jet" cases refer to radio-quiet quasars and 
radio-loud quasars/radio galaxies, respectively. For the latter cat- 
egory, it is unlikely that much of the power in the jet will be de- 
posited inside the galaxy, as essentially all the energy is transported 
out into the lobes which can be located anything from tens of kilo- 
parsecs to a few Megaparsecs away from the central galaxy. The 
additional efficiency factor efr-h accounts for this effect. 

It is not yet certain if and how radio-quiet quasars actually 
provide energy to heat the ambient medium. One possibility is that 
a certain fraction of the b olometric luminosi ty couples to the gas 
through line driving (e.g.. ' Fabian et alj|200 ^. Another possibility 
is that the quasar launches a wind that heats the gas through its ki- 
netic power. It is likely that this process is also linearly coupled to 
the bolometric luminosity. According to the model used in this pa- 
per, low accretion-rate objects are inefficiently accreting up to the 
critical accretion rate. Accretion above that rate is radiatively ef- 
ficient. The approximate dependence of the bolometric luminosity 
on accretion rate is given in eq. il4\ . It has been argued that roughly 
5% of th e bolometric luminosit y is available to heat the surrounding 
gas (e.g.. lDi Matteo et al.l2005l) . We will use this model to compare 
the effects of radiative feedback to the kinetic heating due to jets, 
i.e., the radiative feedback is given by 

-Pradfocd = 0.051/bol- (16) 

Figure [TT| shows the power available from a single source 
(eqns.[75]and|I6l) via both feedback mechanisms, radiative and ki- 
netic (setting epR-ii — 1 for jet-active radio-loud quasars/radio 
galaxies), as function of accretion rate. At low accretion rates, the 
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Figure 11. Power available for feedback as function of accretion rate, both 
in Eddington units (from eqns. |15| and |16) . Sources are kinetically domi- 
nated at low accretion rates, but radiative feedback may dominate at high 
accretion rates. A simple linear dependency may correctly describe the total 
available power as function accretion rate within measurement uncertain- 
ties. 



kinetic feedback of the jet dominates the total power output. At high 
accretion rates, the radiative feedback is efficient and contributes 
significantly to the total power. The jet is either quenched and the 
power in the jet is negligible compared to the power available via 
radiative feedback - or the source is still jet-active and radiative 
and kinetic feedback are comparable for our parameters. If most of 
the energy in the jet is injected into distant lobes and not deposited 
inside the galaxy and epR-n ^ 1, as is very likely the case, the 
kinetic feedback of radio-loud quasars is actually much lower than 
the maximum available feedback energy shown in our figures. 

At low accretion rates, the total power available for feedback 
is linearly dependent on the accretion rate and dominated by the 
kinetic power of the jets. At high accretion rates, the total power is 
either dominated by radiative feedback, or a combination of radia- 
tive and kinetic power. Again, the total feedback depends linearly 
on the accretion rate with a proportionality factor similar to that in 
the hard state. Thus, the assumption that the energy available for 
feedback is roughly Pfccdback ~ O.OlAfc^, seems to be correct 
within a factor 2 for all accretion rates; i.e., Pfccdback is 10% of 
0.1Mc^the power liberated by the accretion process. 



For both AGN and X-ray binaries, both the accretion rate and the 
jet power can be estim ated from either the core radio luminosity 
jKording et"^l2006bl) or the extended low-frequency radio lumi- 
nosity l lWillott et al Il999^. Using these accretion rate estimates, we 
construct accretion rate functions (ARFs) of jet-emitting sources. 
The luminosity function (LF) of low-luminosity sources - which 
likely have slow jets - is roughly in agreement with the bolomet- 
ric luminosity functi on (BLF) of AGN (Fig. ^ as determined by 
iHopkins et al. U2007l) . The ARF of RLQ-like jet sources is roughly 
1 dex below the BLF (Fig. [7}. This provides a measurement of the 
radio-loud fraction with a direct physical meaning: the ratio of vol- 
ume densities of radio-loud and radio-quiet sources at a given ac- 
cretion rate. However, the exact value is unfortunately strongly de- 
pendent on the normalizations used to derive accretion rates. 

We have developed a simple model based on the universality 
of accretion physics in XRBs and AGN that reproduces both our 
ARF as well as the BLF. At low luminosities, all sources are radia- 
tively inefficient and the luminosity depends quadratically on the 
accretion rate (eg. 1141). Sources are radiatively efficient only at high 
luminosities (> 3% Eddington). In our model, we assumed that 
the distribution of accretion rates does not depend on the black- 
hole mass and can be written as a simple broken power law for the 
Eddington-scaled accretion rate (eq. |I3b . With this model we can 
reproduce the measured ARF as well as the local BLF (Fig. llOb. 
Thus, this model can solve the problem that there seem to be too 
few LLAGN compared to the number of weak quasars. In this sim- 
ple model the break in the luminosity function is due to the ineffi- 
cient accretion at low accretion rates (compare Jester 2005). 

Using the corresponding jet power measures, we have calcu- 
lated kinetic luminosity functions and their cosmological evolution. 
Our findings support the idea that the majority of the kinetic power 
is created by lower-luminosity AGN with their likely slow jets. The 
highly luminous radio-loud quasars do not contribute significantly 
due to their low number density. The total power injected into the 
ambient medium by jets from AGN is comparable or even exceeds 
the effect of supernovae (the situation may and will be different for 
individual objects). Only at high redshifts, "quasar-mode" feedback 
mechanisms provide a comparable amount of energy for heating 
the interstellar/intergalactic gas. Finally, we discuss the effects of 
the different feedback mechanisms and find that a roughly constant 
fraction (5-10%) of the accretion power is available for feedback, 
independent of the nature of the AGN. 
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